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Integrated optics allow the generation and control of increasingly complex photonic states on chip 
based architectures. Here, we implement two entangled qutrits - a 9-dimensional quantum system 
- and demonstrate an exceptionally high degree of experimental control. The approach which is 
conceptually different to common bulk optical implementations is heavily based on methods of 
integrated in-fiber and on-chip technologies and further motivated by methods commonly used in 
today’s telecommunication industry. The system is composed of an in-fiber source creating entangled 
qutrit states of any amplitude and phase and an on-chip integrated general Multiport enabling the 
realization of any desired local unitary transformation within the two qutrit 9-dimensional Hilbert 
space. The complete design is readily extendible towards higher-dimensions with moderate increase 
in complexity. Ultimately, our scheme allows for complete on-chip integration. We demonstrate the 
flexibility and generality of our system by realizing a complete characterization of the two qutrit 
space of higher-order Einstein-Podolsky-Rosen correlations. 


INTRODUCTION 

Higher-dimensional entangled quantum states are 
highly relevant for fundamental questions in quan¬ 
tum physics and have gained increasing practical rele¬ 
vance due to their distinguishing properties compared 
to qubit states [1-5]. Moreover, there is a wide range 
of open theoretical questions for quantum systems not 
composed of multiple qubits [6]. The reasons for the 
difficulty in implementing higher-dimensional systems 
are mainly technical in nature. In particular, when 
path encoded photonic states are used phase stabil¬ 
ity quickly becomes very challenging using a bulk op¬ 
tical approach, traditionally used in many quantum 
optical experiments. Consequently, bulk optical de¬ 
signs are restricted in their complexity thereby hinder¬ 
ing progress towards higher-dimensional quantum sys¬ 
tems. Here we choose a conceptually different exper¬ 
imental approach. We utilize recent progress within 
the area of optical telecommunication in management 
and automatization of high numbers of complex op¬ 
tical devices. In particular, we aim at combining 
state-of-the art optical fiber network and integrated 
photonic circuit technology with methods of quantum 
optics. The latter has already proven successful in 
different quantum experiments utilizing various inte¬ 
grated photonic devices [7-13]. Eventually, an inte¬ 
grated optics approach allows for complete access to 
complex quantum systems with a high degree of flex¬ 
ibility and automatization, good scalability towards 
higher-dimensional realizations, as well as direct con¬ 
trol over the parameters defining the quantum system. 
Reconfigurable integrated quantum circuits for qubits 
have been shown so far in a silica-on-silicon architec¬ 


ture [10]. We are not only going significantly beyond 
this work by demonstrating universal unitary opera¬ 
tions in two local 9-dimensional Hilbert spaces, but 
thereby also show the scalabilty of our approach to¬ 
wards even higher-dimensional states. Moreover, we 
achieve interfacing our universal multi-port chips with 
a maximally path-entangled qutrit source - something 
that we demonstrate here for the first time and show 
how the source-design, phase-stabilization and inter¬ 
facing scales very favorable for extending it to higher 
dimensions. Lastly, we enter a new wavelength regime 
for complex integrated quantum circuits, by operating 
in the telecommunication band around 1550 nm. This 
makes our approach directly compatible to fiber-optic 
networks and quantum communication applications 
and is technologically highly promising, as this wave¬ 
length regime features extremely mature fiber tech¬ 
nology and is the main focus of technological advance¬ 
ments in integrated optics. In the context of higher¬ 
dimensional quantum systems, higher-order non-local 
correlations between two (entangled) particles are of 
importance as they represent a fundamental distin¬ 
guishing property to the classical world [14]. However, 
controlled creation and transformation in addition to 
N outcome projective measurements of the higher¬ 
dimensional entangled system are required. All these 
requirements are met by our implemented system. We 
will present the concepts of our approach along with 
a first realization of a fully controlled and automated 
entangled two qutrit 9-dimensional quantum system. 
The final setup allows us to completely map the two- 
qutrit correlation space and to compare it to theoretic 
predictions. 
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(D) Coincidence detection 


FIG. 1: Concept of an entangled two quNit N 2 - 
dimensional quantum experiment. A source (S) creates 
the path-encoded entangled two quNits state. After sepa¬ 
ration each quNit passes N — 1 phase shifters followed by 
a local unitary transformation (U). Eventually the state 
is detected and coincidence detection events are recorded 
(D). 


correlations coincide with the number of all theoret¬ 
ically possible correlations [15]. Naturally, the ques¬ 
tion arises of what is to come for N >2. In contrast 
to a qubit system, a single qutrit (TV = 3) cannot 
be explained by non-contextual approaches [17]. For 
N > 4 infinite classes of symmetric beam splitters can 
exist with consequences for the observation of higher- 
order entanglement [15]. At N = 6 the question of 
the existence of a complete set of mutually unbiased 
basis is still unanswered [18]. Our motivation is to 
obtain experimental access to such quantum systems 
providing the basis for future tests and verifications 
of theoretical concepts. 

It is worthwhile mentioning a particular experiment 
can be considered as one element within the full range 
of possible experiments covered by the general con¬ 
figuration in Figure 1. In contrast, we do not aim 
towards a single predefined measurement but rather 
towards a general platform allowing for a great range 
and flexibility within the space of all possible experi¬ 
ments covered by the configuration in Figure 1. 


A general quantum-optical configuration is illus¬ 
trated in Figure 1. For N = 2 it represents the classic 
text book case of a basic two qubit (N = 2) quantum 
system. Generally the system is composed of three 
components: the creation (S), transformation (U) and 
detection (D) of a quantum state. A source (S) creates 
two particles each composed of N modes, a so-called 
quNit, forming the entangled state 

N 

\iP N ) = (1) 

fc=l 

After separation, each quNit A and B is transformed 
by a local unitary transformation (U): 

W n )->(U a ®UbMn)- ( 2 ) 

The experiment concludes with the measurement e.g. 
single photon detection of both quNits A and B cor¬ 
responding to an underlying projective measurement 

Pab = |a)(o| ® \b)(b \, (3) 

of the quantum state | ip) . The expectation value 
(P a b) relates to the number of detected photon pairs 
and represents the measurement result. 

Although, for N = 2 Figure 1 becomes a conceptu¬ 
ally simple example it has revealed a number of quan¬ 
tum effects with consequences ranging from theoret¬ 
ical insights to prospects of applications. Quantum 

entanglement, manifesting itself in so-called Einstein- 
Podolsky-Rosen (EPR) correlations, is the first ex¬ 
ample to name [15]. Interestingly, although different 
sets of EPR correlations exist in any dimension N, 
only for IV = 2 does the number of different perfect 


EXPERIMENTAL IMPLEMENTATION 

In principle, a photonic quNit system can be en¬ 
coded using any degree of freedom. However, physi¬ 
cal and technical constraints pose a limit. Here, path¬ 
encoding will be utilized. In some way it may be con¬ 
sidered as the most direct way of encoding quantum 
information into a system. Of immediate advantage 
is the possibility of near perfect (experimental) sep¬ 
arability of path-encoded modes via spatial separa¬ 
tion. Additionally, there is no fundamental limit on 
the number of paths. Choosing path-encoding as the 
starting point methods of on-chip and fiber optics of¬ 
fer exceptional mode and path control. The system 
is built using only integrated components and can be 
deconstructed into three parts: Source (S), transfor¬ 
mation (U) and detection (D) as illustrated in figure 1. 

Source of entangled qutrits (S) 

We have previously presented a scalable in-fiber 
source of path-encoded entangled quNits at telecom 
wavelengths [19]. It uses N non-linear crystal waveg¬ 
uides to create photon pairs of same polarization but 
different wavelength centered around 1550nm. All 
crystals are coherently pumped by a common 775nm 
pump beam, split by a sequence of beam splitters. 
The result is a superposition of pair creation events 
occurring at one of the N crystals. This is followed by 
100GHz in-fiber filters separating and filtering each 
photon pair by wavelength (1551.7nm, 1548.5nm). A 
number of optical in-fiber devices ensure the match- 
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ing of optical properties (spectrum, coherence, ar¬ 
rival time, loss, polarization), making all pair creation 
events indistinguishable from each other. Thereby an 
entangled two quNit state is obtained. The main fea¬ 
tures of this design are its high brightness and quality 
of the entangled state, high robustness due to a fully 
fiber-integrated approach, compatibility to standard 
path-encoded integrated optical circuits along with a 
low propagation loss and good scalability in complex¬ 
ity when increasing the dimension of the system. We 
have extended the source of two entangled qubits fol¬ 
lowing the scheme previously presented [19] to two 
entangled qutrits (N = 3). 

The Multiport (U) 

The so-called general Multiport (MP) allows an ex¬ 
perimenter to realize any unitary transformation de¬ 
fined by TV 2 — 1 parameters. Here it is realized using a 
combination of phase shifters and beam splitters. The 
setting of each optical device directly relates to a pa¬ 
rameter of the unitary transformation. Depending on 
the phase and reflectivity settings any unitary trans¬ 
formation is realized [15], [20]. In order to effectively 
manage the internal complexity of the device we have 
chosen an integrated photonic implementation. 

Light-guiding Silica waveguides of low index con¬ 
trast (A ss 1.5%) optimized for C-band operation 
(cut-off wavelength A cutoff ~ 1450nm) are used. 
Waveguide modes closely match to the mode of a sin¬ 
gle mode fiber at 1550nm. Beam splitters are realized 
by evanescent coupling of two waveguides brought to 
proximity. Two 50/50 beam splitters form a Mach- 
Zehnder interferometer with phase shifters in each 
arm representing a tunable beam splitter of arbitrary 
reflectivity. Phase shifters are realized via the thermo- 
optical effect using heating electrodes above the pho¬ 
tonic circuit on top of the cladding layer. Trenches are 
etched around phase shifters to ensure high thermal 
insulation. All phase shifters are permanently con¬ 
nected using a custom electronic probe card mounted 
from the top onto the chip. The 3x3 MP circuit fea¬ 
tures 6 beam splitters and 21 phase shifters (not all 
theoretically required). The chip is thermally stabi¬ 
lized ensuring optimal operating temperature. Total 
insertion loss including in- and out-coupling amounts 
to approximately 9dB. The length of the 3x3 MP cir¬ 
cuit is 6.5cm at a waveguide spacing of 250pm. The 
device is permanently packaged to standard optical 
fibers and electric connectors serving as a convenient 
optical and electrical interface. The complete housing 
has a size of less than 10cm x 5 cm x 5 cm. Any unitary 
transformation entered by the user is automatically re¬ 
alized by a computer program translating the entered 
unitary matrix down to independent heating currents. 


Switching the Multiport to a unitary transformations 
requires t < 100ms. However, we expect the actual 
switching time to be significantly lower. Compara¬ 
ble Silica circuits based on the thermos-optical effect 
have reached switching times of on the order of 1 ms 
[21], For a single Silica switch switching times on the 
order of 50ps have been reported [22]. Nonetheless, 
the current switching time is negligible compared to 
common single photon detection integration times of 
on the order of 10 seconds. Therefore, switching the 
Multiport does currently not limit the experimental 
procedures in any way. We have thus not further fo¬ 
cused on the issue. Due to the high level of control 
and automation the MP can be considered a black¬ 
box device. Any unitary transformation entered by 
the experimenter is automatically realized. 

Detection (D) 

The outputs a = i of qutrit A are detected by a 
free-running, high performance, low noise single pho¬ 
ton detector. A detection event will simultaneously 
trigger three single photon detectors at output modes 
b = /, 77, III of qutrit B opening the detection win¬ 
dow for 2.5ns. A detection event at B will be treated 
as one detected photon pair leading to the count rate 
CC a b . Output modes (a, b) = (*,/),(*,//),(*,///) 
are measured simultaneously yielding the correspond¬ 
ing probabilities of joint detection P a b oc CC a b- Gen¬ 
erally, CC a b depends on all parameters defining the 
quantum system. In particular, the relative phases 
CCab = CC a b{<p i — i — <j>i) are of importance. 
Note, for a ^ i no physical detectors are necessary as 
explained later. 


PREPARATION AND PERFORMANCE 

The complete setup is illustrated in Figure 2. All 
free parameters defining the quantum system are un¬ 
der experimental control. The Source (S) handles 
the creation of the entangled state: 

• The pump power within each crystal controls all 

amplitudes Ax, A 2 , A 3 . Similar to the splitting 
ratios the coherence of the pump beam could 
be controlled by inserting fiber delay lines for 
other types of quantum states. Here, without 
exception, we set the amplitudes to an equal su¬ 
perposition state Ax = A 2 = A 3 = . 

• An active PID feedback mechanism ensures 
complete control and stabilization of all phases 
0 i, </>!, </> 2 , (f> 2 , 03 , <f >3 of the entangled state (de¬ 
tails in [19]). 
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FIG. 2: Illustration of the experimental arrangement. All components are directly interfaced via standard optical fiber 
connectors or directly spliced together for minimal optical loss. (S) A superposition between N photon pair creation 
events is produced by coherently pumping N non-linear crystals. For each channel the photon pair is separated by 
wavelength followed by a number of in-fiber components ensuring matching optical properties. Eventually, an entangled 
two qutrit state is obtained. Using the interference signal of reference light (A 765 ,A 78 s) inserted backwards into the 
system a PID controller is actively stabilizing all phases of the in-fiber system. (U) Each qutrit is directly connected to 
the packaged Multiports performing any arbitrary local unitary transformation onto the system. (D) Eventually, using 
single photon detectors the measurement concludes when pair detection events are recorded. 


The final state should approximately assume the form 

N 

\^ N )=Y J A k e^~^\kk). (4) 

fc= i 

Two Multiport (U) devices are connected to 
qutrit A and B realizing a local transformations taking 
the form: 

((a/b) ii (a/b ) 12 (a/b) 13 \ 

Ua/b = (a/b) 12 (a/b) 22 e 4 ’ 22 (a/b) 2 z.e 4 ’™ , (5) 

\(a/b) i3 (a/b) 23 e ~' i>23 ( a/b) 33 e ^ 33 ) 

with a total of 16 real parameters. Each Multiport 
can be independently configured to realize different 
local unitary transformations within SU(3). 

Single photon detection (D) completes the exper¬ 
iment. A unitary transformation Up ro j of the detec¬ 
tion basis allows for the realization of any projective 
measurement. Thus, if a certain projective measure¬ 
ment is desired, the Multiport is set to include the 
appropriate transformation Up ro j on top of the trans¬ 
formation of the particular experiment U -A UUp ro j . 
The existence of noise in the system is dominantly a 
result of accidental coincidence detection events due to 
a finite detection time window, where the coincidental 
clicks do not originate from a produced photon pair. 
Noise of this type can be quantified and is subtracted 
for all following measurements and quantities if not 


stated otherwise. Due to the property of the noise, it 
can be reduced by decreasing the created pair rates 
per crystal. 

In anticipation of future experiments, instead of 
measuring a particular fixed set of relative phases 
CC a b(cj >i — i — </>[) we choose a different ap¬ 

proach: all relative phases are scanned stepwise within 
0, • • • ,27r . Thereby, a complete map of the correla¬ 
tion space CC a b is obtained. Afterwards, any par¬ 
ticular set of phases can be extracted based on a fit 
obtained from the complete scan of the space. Any 
measurement herein will utilize this approach. 

Before continuing with more complex measurements 
in the full two-qutrit quantum space and for a better 
understanding of the system, we present some more 
straight-forward measurement results: 

• The previous two-qubit source achieved a state 
fidelity of F = (99.44 ± 0.06)% [19]. The cur¬ 
rent qutrit system involves three two-qubit sub¬ 
spaces. Using the Multiport all subspaces A,B,C 
can be accessed and characterized in a simi¬ 
lar way (see [19]). Compared to the theoreti¬ 
cally expected perfect state we reach a state fi¬ 
delity of F a = (96.22 ± 0.73)%, F B = (93.33 ± 
0.88)%, F c = (96.86 ±0.69)% 

Naturally, the resulting fidelity is lower than the 
optimized qubit source. The qutrit system has 
not been optimized for measurement in a qubit 
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subspace. In fact, the remaining part of the 
unitary space the MP cannot realize are per¬ 
fect qubit transformations excluding one mode 
from any interaction. Consequently, any mea¬ 
surement using such transformations will be of 
less quality when assuming a perfect transfor¬ 
mation. We have not further addressed this is¬ 
sue since qubit subspaces will be of little interest 
building a qutrit system. In any case, a sufficient 
quality of the state is achieved. 

• Extrapolating the two qubit subspace measure¬ 
ments to the full Hilbert space, one immediately 
reveals the potential of the Multiport to sustain 
and manipulate genuine higher dimensional en¬ 
tanglement. With the achieved subspace fideli¬ 
ties and an equal population of every path we 
can use the results of Ref. [ 6 ] in order to show 
that the Multiport is capable of maintaining a 
genuine qutrit entangled state. We find a corre¬ 
sponding lower bound on the number of ” e-bits”, 
i.e. entangled bits, the universal currency of bi¬ 
partite entanglement of E 0 p{p ) > 1.15 which 
is significantly above the maximum obtainable 
value of an entangled qubit system, as qubits 
can never encode more than one e-bit. Further¬ 
more reaching the upper bound of qubits would 
require perfectly prepared pure states and noise¬ 
less channels and measurements. 


exit mode. 

For a detailed analysis, we measure a complete map of 
the two qutrit correlation space CC a b{<t>xi fiy) spanned 
by the two relative phases. Iteratively, <p x is set to 
0°,10°,20°,- • ■ ,360° in steps of 10°. For each value of 
<f> x , c/)y is measured in 30 steps between 0 and 27 t. A 
fit f{(j) y ) is applied to the recorded correlation sig¬ 
nal CCab ^(o ,27r ) w ] rl ji 0 (j) x j s kept constant. Even¬ 
tually, the resulting 36 fitted signals f((j) y ) are dis¬ 
played for each 4> x and for every detector combina¬ 
tion ( i,II ), (*,///)} as slice by slice in Figure 

3. A total of 1080 different measurement points with 
8 seconds integration time have been recorded. At 
maximum a coincidence count rate of around 600/8s 
at a SNR > 10 is achieved. The experiment was 
calibrated, controlled, stabilized, and performed fully 
automatically. The measurement results (Figure 3) 
reveal excellent agreement with theoretic predictions. 
Notably, given the experimental resolution even the 
minima are clearly visible. Some particular observa¬ 
tions are listed below: 

• The correlation pattern (Figure 3) exhibits three 
distinct correlation maxima at the predicted 
phase setting (colored dots in Figure 3) [15]. For 
example, if detector combination («, I) exhibits 
a maximum the other two combinations (i,/J) 
and ( i,III ) show a minimum in detected coin¬ 
cidences. 


RESULTS 


For a first demonstration, we perform a complete 
automated characterization of the higher-order two- 
qutrit EPR correlations [15]. 

For its observation a maximally entangled two qutrit 
state is required along with two 3-dimensional beam 
splitters and local 3-outcome measurements pointing 
at the difficulty in observing higher-order correlations. 
Both Multiports are used to realize the extension of 
the balanced beam splitter to higher dimensions, the 
so-called Bcllport, taking on the form [15]: 





G *) 


( 6 ) 


Note that the matrix elements of the Bellport may 
be viewed as elements of the discrete Fourier transfor¬ 
mation (FT) of the inputs. Thus, if a particle enters 
one input it will exit the device at any unpredictable 
output with equal probability (in analogy to the FT 
of a delta peak). However, in contrast to classical 
physics, quantum correlations of two entangled parti¬ 
cles each entering the similar device allow to impose 
certain constraints and predictions onto the particle’s 


• In principle, the number of possible correlations 
between detector combinations equals Nl = 6 
[15]. However, given the experimental resolution 
(Figure 3) we observe three correlation max¬ 
ima at three distinct phase settings equal to 
the system’s dimension N = 3. The remaining 
N\ — N = 3 possible correlations between detec¬ 
tors do not appear, as is expected from quantum 
theory [15]. 

• Examining the middle graph of Figure 3, the 
maximum correlation appears at (60°, 120°) 
while minima appear at (180°, 360°) and 
(300°, 240°) (as illustrated by the arrows). 
The relative distance is given by (0°,0°) , 
(240°, 120°) and (120°, 240°), as expected. 

• For a qubit beam splitter, a measurement of the 
second output is equivalent to a measurement 
of the first output plus an input phase offset 
of 180° [15], [23]. Consequently, measuring a 
certain point (<j> x ,<f> y ) for all possible 9 detector 
combinations should be equivalent to measuring 
the three points (</> x , (j) y ), (4> x + 120°, <f> y + 240°) 
and {cj) x + 240°, (j) y + 120°) while recording only 
joint detection events between (i,I), (i,//) and 
(■*,///). Figure 3 strengthens this assumption 
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FIG. 3: A full scan of the two qutrit correlation space spanned by the two relative phases CC a b(<l>x, (j>y) is performed 
(top row) and compared to theoretic predictions (bottom row). Iteratively, (f> x is set to 0°,10 o ,20°,- • • ,360° in steps of 
10°. For each value of <j) x , 4> y is measured in 30 steps between 0 and 2 -7T. Eventually, a fit is applied to the recorded 

correlation signal CC a b\ r ^ e< '°’ 2n ' > ■ The resulting 36 fitted signals f{(j>y ) are displayed for each <j> x and for every detector 
combination {(*, I), (i, II), ( i , III)} as slice by slice. The maximum coincidence count rates reach approximately 600/8s 
at a SNR > 10. The overlays of the top right picture are used to highlight special characteristics of the correlation space 
as described in the text. The experiment is calibrated, controlled, stabilized, and performed fully automatically. 


which is the reason for the implemented de¬ 
tection scheme using only a single detector for 
qutrit A. 

• The observed correlation pattern is similar to 
the expected classical interference pattern of a 
three-path Mach-Zehnder-interferometer as pre¬ 
viously observed [24], Here, the pattern is non- 
classical: it becomes visible solely in coincidence 
detection. Local observation of single detection 
events (ideally) does not exhibit any pattern and 
a photon will be observed with the same proba¬ 
bility of 1/3 at any output. 

• Three path interference results in a higher phase 
sensitivity S oc \dl/dcp\ due to the existence of 
side lobes [24]. In order to reveal the effect 
within the two particle correlation pattern the 
output intensity / is replaced by the number of 
detected pairs and extracted from the measure¬ 
ment data. The maximal measured phase sensi¬ 
tivity averaged over the three measured detector 


combinations becomes: 

S = (0.703 ± O.Ol^rad -1 (7) 

The theoretical maximum sensitivity of three- 
paths interference is given by S^ h « 0.782 . 
Compared to the qubit case with S ^ = \ an in¬ 
crease in sensitivity by a factor of 7 = 1.41±0.04 
is achieved. 

These observations demonstrate the good operation 
of the Multiport devices and the achievable purity of 
the input states. The fact that the whole correlation 
space is accessible also implies that any quantum in¬ 
formation protocol, that can be performed using en¬ 
tanglement and local unitary transformations, is di¬ 
rectly implementable with our setup. 

In addition, we wish to highlight the scalability of 
the current system towards higher-dimensional quan¬ 
tum systems. The scalability of the source has been 
previously addressed [ .9]. While a detailed discussion 
would be beyond scope of this article, we can provide 
a representative example of the power of integrated 
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optics in the context of the Multiport design. The in¬ 
tegrated Matrix Switch, a device used in the area of 
communication, redirects light entering one of N in¬ 
puts to one of N outputs by a combination of tunable 
MZIs. Matrix Switches up to 32 x 32 in/outputs based 
on similar Silica waveguide technology have been re¬ 
ported [25] featuring 2048 beam splitters, 2048 phase 
shifters and 60cm of total waveguide length integrated 
on a single chip at an insertion loss below lOdB. In 
comparison, the Matrix Switch circuit structure is 
highly similar to a Multiport. In fact, in terms of de¬ 
vice complexity a 32 x 32 Multiport would require half 
the number of optical devices compared to the above 
Matrix Switch 1 . An important aspect of interfacing a 
high-dimensional path-entangled source as ours with 
multi-port chips is the scaling of the resources for 
active phase-stabilization. The method we have de¬ 
veloped for the N=3 case uses back-illuminating the 
complete beam path with two lasers of different wave¬ 
length injected into the output. The interference sig¬ 
nal at 2 of the inputs for both lasers (4 detectors in 
total, see Fig. 1) is then measured and used for sta¬ 
bilizing the phases of the fibers connecting the source 
and the multi-port chip. A careful analysis shows that 
scaling up this design would require only 2 more de¬ 
tectors (one for each wavelength) per additional di¬ 
mension, but no additional lasers. Thus our active 
phase-stabilization scheme scales very favourable and 
contributes to the excellent scalability of our whole 
approach to higher dimensions. Additionally, losses 
of lOdB are handled by the current source and detec¬ 
tion system. Although of relevance for the circuits’s 
design, due to the nature of integrated circuits coher¬ 
ence length and path difference are assumed not to 
be of major issue. However, a major stepping-stone 
towards high-dimensional Multiports may be the de¬ 
creasing process control with increasing circuit sizes 
and therefore reduced quality of the optical devices. 
Here, the major concern is the extinction ratio or vis¬ 
ibility of the MZIs of the MP. Ultimately, low extinc¬ 
tion ratios will reduce the MP’s capability to realize 
arbitrary unitary transformations. In the following 
we provide an estimate of this effect: The range of 
unitary operations implementable by the MP can be 
quantified. We define the ’’coverage” as the amount 


1 Unfortunately, current-generation Matrix Switches are highly 
optimized leading to a circuit design not usable as a Multi- 
port. Other designs [27] have a higher similarity to the Mul¬ 
tiport, we will address this issue in a separate article. 



FIG. 4: Estimated coverage of the unitary space of an 
N x N Multiport of standard design [1, 15]. Each plot¬ 
ted line assumes a different extinction ratio for all MZI 
units inside the Multiport. Phase shifters are assumed to 
cover the full range of phases as has been experimentally 
verified. Currently, the characterized Multiport exhibits 
an average visibility of (96.3 ± 0.01) % (extinction ratio of 
(14.2 ± 0.1)dB). However, the current MP implementation 
has the a first/prototype run. Consequently, we are con¬ 
fident that next-generation devices tailored towards high 
extinction ratios will provide significant improvements. As 
an example, a 10-dimensional Multiport with an extinction 
ratio of —30dB of all MZIs can realize more than 92% of 
the full unitary SU( 10) space. 


of unitary space ’’covered” by the Multiports: 


COVn = 


max Ai 


n i= i i d\iJ N 

min Xi 

fu(N) dU N 


( 8 ) 


i.e. the ratio between the reachable unitary space 
given by all parameters and the full unitary space. 
The factor Jjv represents the absolute value of the 
Jacobi determinant, see [26] for further details. The 
limits of Xi are related to the range of reactivities 
(phases) implementable by each beam splitter (phase 
shifter). Analyzing a single Multiport we estimate the 
coverage to (96.3 ± 0.01) % of all unitary operations. 
Here, the limitation stems from the extinction ratios of 
the Mach-Zehnder interferometers (MZI) while phase 
shifters cover the full range 0...27T. Using the exper¬ 
imental results we can further extrapolate the cover¬ 
age of a higher-dimensional Multiport as illustrated in 
Fig. (4). Our results thus demonstrate that with ex¬ 
perimentally feasible noise levels, moderate scalability 
of the devices is straightforward. Significant improve¬ 
ments can be expected with further tailoring of next 
iterations of Multiport chips. Note, however, using a 
Multiport of high coverage might not even be neces¬ 
sary as long as it can implement the sets of unitaries 


















required for the experiment or quantum information 
protocol. 

CONCLUSION 

Combining previously developed devices we have 
successfully implemented a system allowing for an ex¬ 
ceptionally high degree of experimental control of an 
entangled two qutrit, 9-dimensional photonic quan¬ 
tum system. An in-fiber source creates an entangled 
two qutrit state of any amplitude and phase, while two 
Multiport devices realize any desired arbitrary local 
unitary transformation onto the quantum state. We 
demonstrate the quality of the system by performing a 
complete characterization of high-order perfect EPR 
correlations between two qutrits. Results are com¬ 
pared to theoretic prediction. Based on the perfor¬ 
mance of the current implementation, no direct tech¬ 
nical limitations towards higher-dimensional systems 
become apparent. The source is directly extendible 
to two ququarts and could be implemented on-chip 
as well [19]. Compared to other state-of-the-art inte¬ 
grated circuits used in the area of telecommunication 
[28] higher-dimensional realizations of the Multiport 
are highly promising for quantum systems of dimen¬ 
sion N. Further experiments may involve the genera¬ 
tion and manipulation of entangled states [29], state 
discrimination [30], simulation of quantum logic [31], 
general optical simulation of spin-1 systems [15], an 
optical analog of the Stern-Gerlach experiment [15], 
tests of non-contextuality [4], [15], realization of gen¬ 
eral POVM measurements [30], [32], or assist in the 
search for the existence of complete sets of MUBs [33] , 
[34], Additionally, the compatibility of the in-fiber 
source and Multiport to standard fiber networks at 
the Telecom wavelength will be of further interest for 
future applications. The design is a viable approach 
towards higher-dimensional quantum systems and of¬ 
fers great potential towards a higher level of integra¬ 
tion and complexity of quantum systems. 
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